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Abstract 

We assume that current state of the Universe can be described by the 
Inert Doublet Model, containing two scalar doublets, one of which is re- 
sponsible for EWSB and masses of particles and the second one having no 
couplings to fermions and being responsible for dark matter. We consider 
possible evolutions of the Universe to this state during cooling down of the 
Universe after inflation. We found that in the past Universe could pass 
through phase states having no DM candidate. In the evolution via such 
states in addition to a possible EWSB phase transition (2-nd order) the 
Universe sustained one 1-st order phase transition or two phase transitions 
of the 2-nd order. 



1 Introduction 

According to the standard cosmological model about 25% of the Universe is 
made of Dark Matter (DM). Different candidates for DM particle are now dis- 
cussed in the literature. One of the widely discussed models is the Inert Doublet 
Model (IDM) [T] - a Z2 symmetric 2HDM with a suitable set of parameters. 
The model contains one "standard" scalar (Higgs) doublet (j>s, responsible for 
electroweak symmetry breaking and masses of fermions and gauge bosons as in 
the Standard Model (SM), and one scalar doublet, <f>n, which doesn't receive 
vacuum expectation value (v.e.v.) and doesn't couple to fermionfl 

In this model four degrees of freedom of the Higgs doublet 4>s are as in the 
SM: three Goldstone modes become longitudinal components of the EW gauge 
bosons and one mode becomes the Higgs boson (here denoted as hs)- All the 
components of the scalar doublet cj)D are realized as massive scalar D-particles: 
two charged and two neutral ones Dh and Da- By construction, they 



1 Our notations are similar to those in the general 2HDM with the change (pi — > cj>g, 02 
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possess a conserved multiplicative quantum number and therefore the lightest 
particle among them can be considered as a candidate for DM particle. 

Assuming, as usual, that DM particles are neutral, we consider such variant 
of IDM, in which masses of £>-particles are 

M d ±,M Da M Dh or M d ±,M Dh ^ M Da . (1) 

Possible masses of these D-particles are constrained by the present accelerator 
and astrophysical data (see e.g. [HG3). 

In this paper we assume that the current state of the Universe is described by 
IDM. We discuss possible variants of the history of the phase states of Universe 
during its cooling down after inflation. In some respects, this analysis can be 
considered as particular case of analysis [J] , [5] . We use below some results and 
notations from [4]- [7]. 



2 The Lagrangian 

In this paper we consider an electroweak symmetry breaking (EWSB) via the 
Brout-Englert-Higgs-Kibble (BEHK) mechanism, as described by the Lagrangian 

C = L s g f+L H +C Y U>tAs), C H =T-V. (2) 

Here, Lgj 1 describes the SU{2) x U(l) Standard Model interaction of gauge 
bosons and fermions, which is independent on the realization of the BEHK 
mechanism. In the considered case the Higgs scalar Lagrangian Ch contains 
the standard kinetic term T and the potential V with two scalar doublets <j)g 
and tf>D- The Cy describes the Yukawa interaction of fermions tpf with only one 
scalar doublet <j>s, having the same form as in the SM with the change <f> — > 4>s. 



Potential. The potential must be Z2 symmetric in order to describe IDM. 
Without loss of generality it can be written in the following form 



m ii Os</>s ) +m| 2 W D 4>d ) 



Ai(4^s) 2 + A 2 (^^)2 + A 3 (4&0(4&>)+ 



with all parameters real and with additional conditional 



(3a) 



A 5 < . (3b) 

2 ln the general Z2 symmetric potential the last term has a form ^As^jjc^o^ + Ag^^^s) 2 ! . 

The physical content of theory cannot be changed by the global phase rotation tj> a — > (f> a e tOLa 
(a = S, D). Starting with an arbitrary complex A5 = |As|e lp we select a$ — old = p/2 + 7r/2, to 
get <[3j> with negative A5 = — |As|. 
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The IDM is realized in some regions of parameters of this potential. To study 
thermal evolution, we will consider also other possible vacuum states of such 
potential, at another values of parameters. 

To make some equations shorter, we use the following abbreviations: 

A 345 = A 3 + A 4 + A 5 , fl=-^=. (4) 

v A1A2 

Discrete symmetries. This potential Q is invariant under two discrete sym- 
metry transformations of a Z2 type: 

S: ^4-0 s , fe^fo, SM^SM, (5) 
D: 4> s ^4>s, 4>d^-4>d, SM-^SM, (6) 

where SM denote the SM fermions and gauge bosons. 

We call these transformations S -transformation and D -transformation, re- 
spectively. In the case when vacuum has vanishing vacuum expectation values 
(0s) = {4>d) = the mentioned above invariance of V results in the Z?-parity 
and S'-parity conservation in the processes involving only scalars (or scalars and 
gauge bosons). The Yukawa term violates S'-symmetry even if (<j>s) = (0X>) = 0; 
while it respects D-symmetry in any order of perturbation theory. 

Positivity constraints. To have a stable vacuum, the potential must be 
positive at large quasi-classical values of fields \<pi\ (positivity constraints), for 
an arbitrary direction in the (</>s, </>r>) plane. These conditions limit possible 
values of Aj (see e.g. [8]). In terms of parameters Q positivity constraints 
which are needed in our analysis, can be written as 

Ai>0, A 2 >0, R + 1>0. (7) 



3 Thermal evolution 

Main goal of this paper is to consider an evolution of the Universe during its 
cooling down to the present inert phase. For this purpose we consider thermal 
evolution of the Lagrangian, following the approach presented in [SJ [5] . 

Potential. Since the Hubble constant is small, we assume a statistical equilib- 
rium at every temperature T. In this approximation, at the finite temperature, 
the ground state of system is given by a minimum of the Gibbs potential 

V G = Tr (Ve~ k / T ) /Tr (e~"' T ) . (8) 

In the first nontrivial approximation and high enough temperature the obtained 
Gibbs potential has the same form as the basic potential V§5§, i.e. as the 
potential at zero temperature. The coefficients A's of the quartic terms in the 
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potential Vq and V coincide, while the mass terms vary with temperature T, 
as follows 



m 2 n (T) = mj, - Cl T 2 , m 2 2 (T) = m| 2 - c 2 T 2 , 

_ 3Ai + 2A 3 + A 4 3 g 2 + , 9 ' 2 gj + gf 
Cl 12 + 32 + 8 ' W 

3A 2 + 2A 3 + A 4 3g 2 + <?' 2 

C2 = — 



12 32 

Here 5 and 5' are the EW gauge couplings, gt ~ 1 and ~ 0.03 are values of 
the SM Yukawa couplings for t and b quarks, respectively. 

Generally each of coefficients c\ and c 2 can be either positive or negative. 
However, in virtue of positivity conditions (J7J their sum is positive, 

C 2 + C! > 0, (10) 

even neglecting (positive) contributions from gauge bosons W/Z and fermions. 

We will show later on that for a realization of the present inert vacuum with 
neutral dark matter particle one needs A4 + A5 < (|2"4")l . Therefore, at R > 
we have A3 > 0. Taking into account that A5 < (|3b|) . we obtain that c\ > 0, 
C2 > 0. At R < there are no constraints on signs of c\$'- 

R > : ci > 0, C2 > 0; R < : arbitrary signs of ci,2- (11) 

Yukawa interaction. The form of Yukawa interaction and values of Yukawa 
couplings don't vary during thermal evolution. 



4 Extrema of the potential 

Following [4] we first consider extrema of the potential Q at arbitrary values 
of parameters. The extrema conditions: 

dV/84H\^=^)=0, £>Wlk =(0i> = 0, {i = S,D) (12) 

define the extremum values (4>s) an d {4>d) of the fields 4>s and (f>£>, respectively. 
The extremum with the lowest energy {the global minimum of the •potential) re- 
alizes the vacuum state of the system. Other extrema are saddle points, maxima 
or local minima of the potential. 

The most general solution of (|12p can be written in a following form: 

to>=^(°). <^ = 7iU)' (« 9 =«a+i«2>i+« a ) (13) 

since for each electroweak symmetry violating extremum (EWv) with (4>s) 7^ 0, 
one can choose the z axis in the weak isospin space so that ((f>s) ~ ^ , with 
real, nonnegative vs (choosing a "neutral direction" in the weak isospin space). 
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Neutral extrema. The solutions of (|12p with u = are called neutral ex- 
trema, as they respect U{\) symmetry of electromagnetism. For these extrema 
the conditions ([12")) can be written as a system of two degenerate cubic equations: 

us(-m? 1 + Aiu| + A 345 uf,) = 0, u D (-m| 2 + A 2 uf, + A 345 ?;|) = , 

(14) 

^1 > , v 2 D ^Q. 

This system has four solutions, one solution defines electroweak symmetric 
extremum EWs and three solutions define EWSB extrema: inert extremum I\ , 
inert- like extremum 7 2 and mixed extremum M. Below we list their v.e.v.'s and 
extrema energies £ a : 

EWs : v D = 0, v s = 0, 

h : v D = 0, v 2 s = v 2 = 

1% : v s = 0, v 2 D = v 2 = 



M 



2 _ mjjAg - \345mi2 
MM - A 345 

mf 1 A 2 - 2 

iM = 





^EWs — 0; 


(15) 


Ai ' 


p m ii. 
£h 8A? 


(16) 


A 2 ' 




(17) 


v 2 D = 




m| 2 Aa - A 3 4 5 m? 1 _ 




AiA 2 - A§ 45 
9 4 \ 


(18) 



8(AiA 2 - A 345 ) 



Some of the equations (|T6l) - (fT8|) can give also negative values of v s or vf^, in 
contradiction with basic condition for the extremum (|14j) . In such case the 
extremum, described by corresponding equations, is absent. 

The energy differences between 7i j2 and M extrema are as follows: 

c c (mf 1 A 3 45 - m 2 ^) 2 (m| 2 A 34 5 -m?iA 2 ) 2 

£ ^- £ ^- 8X 2 X 2 (l~R 2 ) ' E ^- Em - 8AiA|(1 - R 2 ) ' {W) 



Charge breaking extremum. For ii^O the extremum violates not only EW 
symmetry but also the U(l) electromagnetic symmetry, leading to the electric 
charge non-conservation. According to general analysis in jTOl [TTJ [T^l 01 [7] this 
extremum can realize vacuum state only if: 

A 4 + A 5 > 0. (20) 

We will see later on that at this condition the DM particle can not be neutral, 
that contradicts ([T]). 



5 Vacuum states 

Below we describe briefly properties of neutral extrema, listed in the previous 
section, provided that they are realized as true vacua. 
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5.1 Electroweak symmetric vacuum EWs 



The electroweak symmetric extremum with (<j>g) = (4>d) = exists for all values 
of parameters of the potential It respects the D and ^-symmetries of the 
potential. This extremum is a minimum, realizing vacuum state, at 

m? x < 0, m 2 2 < 0. (21) 

In this case, gauge bosons and fermions are massless, while scalar doublets 4>S 
and 4>d have masses equal to \m\i\ and 1 17122]) respectively. 



5.2 Inert vacuum l\ 

In the case when I\ extremum realizes vacuum, the Inert Doublet Model de- 
scribes reality. The standard field decomposition near I\ extremum has a form 

I G+ \ ( D+ \ 

(j) S = v + hs + iG , cf> D = D H + iD A , (22) 

V V2 / \ V2 J 

where and G are Goldstone modes, while h$ and D = Dh, Da,D ± are 
scalar particles. Here the Higgs particle h$ interacts with the fermions and 
gauge bosons just as the Higgs boson in the SM. 



Symmetry properties. The inert vacuum state violates the S'-symmetry ([5]). 
However, this state is invariant under the D-transformation ([5]) just as the whole 
basic Lagrangian ([2]). Therefore the D-parity is conserved, and due to this fact 
the lightest Z?-particle is stable, being a good DM candidate. 



Allowed region of parameters. For the inert extremum to exists it is nec- 
essary that m\i > (|16p . In accordance with (fT6|) and (|17[) . the extremum 
l\ can be a vacuum only if m\^/ '\/Ai > m\ 2 / \T^2- Additional condition arises 
from a comparison of I\ and M extrema. In virtue of (|19j) at 1 — R 2 < the 
extremum M can exist but its energy is larger than energy of I\ extremum - so 
that the extremum 1\ realizes vacuum. At 1 — R 2 > the inert extremum still 
can be a vacuum, in the case when the mixed extremum does not exist, i.e. if 
at least one of quantities Vg, v 2 D defined by eq. ([T8|) is negative. 

Note, that due to the positivity constraint 1 + R > in the case when 
1 — R 2 < we have R > 1. For the opposite case, with 1 — R 2 > 0, the quantity 
R can be either positive or negative. 



Particle properties. The quadratic part of the potential written in terms of 
physical fields hs, Dh, Da and D ± (f22|) gives the following masses of scalars: 

M 2 = \iv 2 = m 2 , M 2 



hs - -"'11' 1V1 D± - o ' 

1 (23) 
M 2 =Ml ± + ^—^v\ M 2 =M 2 n± + ^±^v 2 . 
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The requirement that lightest P-particle is a neutral one ((TJ) results in the 
condition 

A 4 + A 5 < 0. (24) 

Since A5 < (|3bf the "scalar" Dh is lighter than " pseudoscalar"0 Da- 

As in the standard 2HDM, scalars Dh and Da have opposite P-parities but 
since they don't couple to fermions, there is no way to assign to them a definite 
value of P-parity. However, their relative parity does matter and for example, 
vertex ZDhDa is allowed while vertices ZDhDh and ZDaDa are forbidden. 
According to our basic assumption on the Yukawa interaction ^ P-particles 
don't interact with fermions. Neither there are interactions of P-particles with 
gauge bosons V of the type PjViV^. 



5.3 Inert-like vacuum I 2 

The inert-like vacuum I2 is "mirror-symmetric" to the inert vacuum p, compare 
(fl6|) and (fl7|) . The interaction among scalars and between scalars and gauge 
bosons are mirror-symmetric as well, so the only difference between p and p 
arises from the Yukawa interaction. 

Main formulae for this state are similar to those for the vacuum p with 
obvious replacements. The corresponding field decomposition is given by 

I S+ \ ( G+ \ 

cj) S = S H + iS A , <Pd=\ v + h D +iG , (25) 

\ V2 J V V2 J 

with one Higgs particle Kb and four S'-particles: Sh , Sa, S ± . 



Symmetry properties. The inert-like vacuum p violates P-symmetry ©. 
This state as well as the Higgs potential is invariant under the 5-transformation 
(|5|). However, in contrast to the inert vacuum, here S'-parity is not conserved 
by the whole Lagrangian because of the form of Yukawa interaction. 



Allowed regions of parameters. The inert-like extremum exists for > 
0. In order to have an inert-like vacuum it is necessary that m\ 1 l\f\\ < 
m 22/ V~^2- For 1 — P 2 < there are no additional demands. If 1 — P 2 > 
inert- like extremum can be a vacuum only if at least one of quantities Vg, ffj, 
defined by eq. (fT8|) . appears to be negative. Both these conditions are similar 
to those for the inert vacuum 1%. 

3 Note that the rephasing transformation (f>± — > <j>i, <j>2 i<t>2, changing sign of A5, results in 
change Dh o Da in I\ state. 
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Particle properties. The masses of the Higgs boson hp and S-scalars are 
given by (cf. (|23]>) 



Jlf2 _ \ ,,2 _ ^2 71^2 „ ^ Jt/ "11 

M hD - A 2 u - m 22 , iW s± - , 

2 (26) 

Ml A = M| ± + , M| fl = MJ ± + . 

The Higgs boson hjj couples to gauge bosons just as the Higgs boson of the 
SM, however it does not couple to fermions at the tree level. The S'-scalars do 
interact with fermions. Therefore, here there are no candidates for dark matter 
particles. That is the reason to call this vacuum inert-like vacuum. 

Note that all fermions, by definition interacting only with 4>s with vanishing 
v.e.v. (4>s) = 0, are massless. (Small mass can appear only as a loop effect.) In 
such vacuum state particles form roughly uniform plasma with massless fermions 
and heavy gauge bosons and scalars. 



5.4 Mixed vacuum M 

The mixed extremun0 M violates both D- and .^-symmetries, i.e. the full Z 2 
symmetry of the potential. In this vacuum we have massive fermions and no 
candidates for DM particle, like in the SM. The decomposition around the mixed 
vacuum looks as follows: 

( pt \ ( pi \ 

4> S = VS + PS + »X5 , <f>D = \ VD+ PD + IXD , (27) 

V V2 J \ V2 J 

where the p~g and p\-, lead to two orthogonal combinations G + and H + , while 
ps and pd (xs and \d) - to two orthogonal combinations h and H (G and 
A), respectively. There are here five Higgs bosons - two charged i? ± and three 
neutral ones: the CP-even h and H and CP-odd A. 



Allowed regions of parameters. In accordance with (fl8|) and (|19l) the 

mixed extremum is global minimum of potential, i.e. vacuum, if and only if 
the following conditions hold: u| > 0, v 2 D > and 1 — R 2 > 0. For v.e.v. 's 
squared given by eqs. (jT5J the latter conditions can be transformed to the rela- 
tions between mass parameters and m 2 2 '- 

at 1 > R > 
at > R > -1 



4 Sometimes called a normal extremum N, see e.g. |12l 



< R 



'22 



> R 




(28) 
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Particle properties. Masses of scalars are as follows (see, e.g. 011]) 

M£ ± =-^V, M 2 a = -v*X 5 , (v*=v 2 s +v%). (29) 
The neutral CP-even mass matrix is equal to 

M=(, Aiw| WsM (30) 

\A345VSVD MVjj J 

Note, that the extremum can be minimum only if both diagonal elements of 
mass matrix and its determinant are positive, i. e. \i\2V^v 2 D (l — R 2 ) > 0, in 
agreement with the above mentioned conditions. It means also that in the case 
if mixed extremum is minimum, it is global minimum - vacuum. 

The mass matrix (|3U|) gives masses of the neutral CP-even Higgs bosons: 

2 Ai^| + X 2 v 2 D ± V(Ai4 + \ 2 v 2 D ) 2 - 4detAl 

m h,H — 2 ' ^ ' 

with sign + for the H and sign — for h. 

Couplings of the physical Higgs bosons to fermions and gauge bosons have 
standard forms as for the 2HDM, with the Model I Yukawa interaction. 



6 Evolution of phase states of the Universe 

In this section we consider possible phase history of the Universe, leading to the 
inert vacuum I\ today, using the thermal evolution described in sec. 3. 

To summarize properties of different vacua of the ^-symmetric potential 
and to classify all possible ways of evolution of the Universe we will use phase 
diagrams in the (pi(T), ^jl 2 {T)) plane, where 

M{T) = m 2 n (T)/^, M2 (T) = m 2 22 {T)/^\ 2 . (32) 

Let us remind (sect. [3]) that in our approximation during cooling down of 
Universe parameters \ are fixed, while mass parameters rn 2 i vary. These varia- 
tions result in modification of vacuum state and a possible change of its nature. 
Possible types of evolution depend on value of parameter R (j4|) and are depicted 
in the figures 16. 11 16.21 and 16.31 The possible current states of Universe are rep- 
resented in these figures by small black dots P = [pti, ^2) 0- Since currently 
we are in the inert phase, we have ^1 > for each point P (see sect. 15. 2| ). The 
parameter /X2 can be both positive (points PI and P3) and negative (points P2, 
P4 and P5). 

In accordance with ^ a particular evolution leading to a given physical 
vacuum state P is represented by a ray, that ends at a point P. Arrows on 

5 ln subsequent discussion we will distinguish present day values of parameters fii = ^;(0) and 
their values fii(T) at some temperature T. 
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these rays are directed towards a growth of time (decreasing of temperature). 
The direction of the ray is determined by parameters (cf. 



ci = c 1 /^/Xl, 



C2 = C2/ \f>V2, 



c = c 2 /ci. 



(33) 



For different possible positions of today's point P we consider typical evo- 
lutions for different possible values of parameter c. In figures below all repre- 
sentative rays are shown; they are labeled by two numbers, with the first one 
corresponding to the label of the final point P. 



Phase diagram for this case is presented in Fig. 16.11 It contains one quadrant 
with EWs phase and two sectors, describing the I\ and 1% phases. These two 
phases I\, I2 are separated by the phase transition line /zi = y,i (thick black line). 
Two typical positions of today's state are represented by points PI (/i2 > 0) 
and P2 (/x 2 < 0). Since (according to (fTTjl ). both c\,C2 > (c> 0), all possible 
phase evolutions are represented by rays 11 and 12 for the today's point PI and 
by a ray 21 which leads to the today's point P2. 



Ray 11: c > /Z2/M1 > 0. The Universe started from the EWs state and after 
the second-order EWSB transition at mf 1 (T) = 0, i.e. at the temperature 



6.1 The case R > 1 




Figure 1: Phase diagram for R > 1 case. 




(34) 



has entered to the present inert phase I\. 
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Ray 12: < c < Hz/fli. The Universe started from the EWs state. Then it 
went through the EWSB second-order phase transition into the inert-like phase 
I2 at m\ 2 {T) — 0, i.e. at the temperature equal to 



T E Ws,2 = Y«^I 2 / C 2 = V^M- (35) 

The next transition is the phase transition from the inert-like phase 
the today's inert phase I\ at the point R, where /^(T") = /xi(T), i.e. 
temperature 

V ci - c 2 

That is the first-order phase transition with the latent heat given by 

= (fi2C 1 -fiic 2 )TlJA. (37) 

^(t) en 

Ray 21: [12 < 0. The Universe started from EWs state and after the second- 
order EWSB transition at the temperature (|34|) has entered the today's phase. 

6.2 The case 1 > R > 



<9£ /2 



dT 



dT 




Figure 2: Phase diagram for 1 > R > case. 

The phase diagram in Fig. 16.21 is modified in comparison with the Fig. 16. 1[ 
as according to (|28|) in the upper right quadrant of the considered (/zi, ^2) plane 
the new sector - the mixed phase M - occurs in the region 

< Rm < H2 < m/R- (38) 

As before, since R > both C1.2 > and consequently c > 0. 
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Since currently we are in the inert vacuum, the possible today's states are 
of type of points P3 and PA, for which 



[l 2 < RfJ-i- 



(39) 



All possible phase evolutions are represented by three rays in Fig. 16.21 with rays 
31 and 32 having the today's endpoint P3 while the ray 41 is pointing PA. 

For the rays 31 and 41, phase evolutions are as for the rays 11 and 21, 
respectively. New situation appears for the ray 32. 

Ray 32: < c < lia/pi. The Universe started from the EWs state. Then at 
the temperature given by (1331) it went through the EWSB second-order phase 
transition into the inert-like phase Ii- At the subsequent cooling down the 
Universe goes through the mixed phase M into the present inert phase I\ . The 
second-order phase transitions I2 — > M and M — > I\ happened at the following 
temperatures 



In accordance with equations in sect. [5j at the transition point 1% —> M masses 
of Sh and h vanish, while at the transition point M ~ > I\ masses of h and Du 
become 0. At small distance from the transition point with temperature T p htr 
these masses grow as a function of the temperature T as M% = A a \T 2 — T^ htr \, 
with different coefficients A a . 



phtr ' 




(40) 



6.3 The case > R > -1 




Figure 3: Phase diagram for > R > — 1 case. 
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The phase diagram is presented in Fig. 16.31 In this follows from 

(128)) . the mixed phase M region is realized in a wider region than in Fig. 16.21 
(for < R < 1), even beyond an upper right quadrant of this plane0, namely: 

M2 > Pi/R, M2 > Mi-R- (41) 

Since currently we are in the inert vacuum, we have 

/i 2 < RfH (P2 < 0). (42) 

Therefore, in this case we have only one type of today's point P5. However, new 
opportunities appear due to larger freedom for temperature coefficients Cj, as in 
accordance with pip in this case C\ and 82 can be either positive or negative. 

All possible phase evolutions leading to the point P5 are represented in 
Fig. 16.31 by four rays 51, 52, 53 and 54. 

The ray 51 describes similar evolution as rays 21 and 41. New are rays 52, 
53 and 54 with common feature, which is a lack of electroweak symmetry in 
very early stages of the Universe. 

Ray 52: c\ > 0, ca < 0, c > /X2/mi- Here a high-temperature state of the 
Universe is the inert-like vacuum 12- With cooling down the Universe goes 
through electroweak symmetric phase EWs into the present I\ phase. The 
second-order phase transitions I2 — > EWs and EWs —t I\ happened, respec- 
tively, at the temperatures 

Ti.ews = \/ M2/C2 , T E v/s,i = x/Mi/ci • (43) 

Ray 53: c\ > 0, C2 < 0, c < ^2/ Hi- Here a high-temperature state of the 
Universe is an inert-like vacuum I-2- With cooling down the Universe passes 
through the mixed phase M into the present 1\ phase. The phase transitions 
I2 —> M and M — > I\ are of the second order; they happened at the temperatures 
given by eqs. (|4T)|) . 

Ray 54: c\ < 0, C2 > 0. For this ray the Universe stays in the inert vacuum 
I\ during the whole evolution. 

6.4 Summary of possible evolutions 

We find that in the considered approximation the thermal evolution of Universe 
to the current inert phase can be studied effectively in the (/ii, ^2) plane, at fixed 
values of quartic parameters \ . Different types of such evolution represented as 
directed rays depend crucially on two parameters: R Q, describing the allowed 
for various vacua regions of the (/ii,/^) plane, and c (|3"3")l . determining the 
direction of rays. The first one depends only on the ratios between coefficients 

e This case was overlooked in the literature, see e.g. [I] and also [2]. We thank G. Gil and 
B. Gorczyca for discussion on this point. 
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of quartic part of potential Aj, while the second depends both on mentioned 
parameters Ai, which are unknown up to now, and on precisely known gauge 
and Yukawa couplings. 

One can distinguish following types of evolution to the current inert phase: 

• I. The Universe evolves from the initial electroweak symmetric state. 

— Ia. The simplest evolution to the inert phase is realized through a sin- 
gle EWSB phase transition of the second-order type (rays 11, 21, 31, 
41, 51). Dark matter appears at this single transition simultaneously 
with EWSB. 

— Ib. After the first EWSB phase transition Universe passes into the 
inert- like phase 1%. Then it passes into the inert phase I\ either 
directly (first-order phase transition, ray 12), or through the mixed 
phase M (two second-order phase transitions, ray 32). In both cases 
dark matter appears only after the last phase transition to the inert 
phase. 

• II. The Universe evolves from the initial state having no electroweak sym- 
metrjQ. 

— Ha. The initial phase is the inert one I\ . Evolution contains no phase 
transitions. Dark matter existed always (ray 54). 

— lib. The initial phase is the inert-like one I^- It contains no dark 
matter. Evolution to the current inert phase 1\ undergoes through 
two second-order phase transitions either via the mixed phase M 
(ray 53), or via the EWs phase, i.e. with a temporary appearance of 
electroweak symmetry (ray 52). In both cases dark matter appears 
only after the last phase transition to the inert phase I\. 

Each of these evolutions can be realized in wide range of parameters. 

7 If DM is charged 

Model independent analysis shows that the case with charged DM particle is not 
ruled out absolutely, but charged DM particles must be heavier than lOOg TeV, 
where q is electric charge of DM particle in units of electron charge [2]. For 
IDM it means M H ± > 100 TeV. Such a heavy mass seems to be unnatural in 
the modern particle physics with natural energy scale < 1 TeV. In this case, our 
new model will include both new particles and new energy scale of phenomena. 
Obviously, such an opportunity cannot be tested at colliders in the estimable 
future. 

7 Such opportunity was discussed by a number of authors — see e.g. | 13 |. Certainly, it is not 
ruled out, but it contradicts a key idea of modern approach - the state at very high temperatures 
has high symmetry which is broken at cooling down of the Universe. In this sense this opportunity 
is unnatural. 
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The case with charged DM particle can be realized in IDM only if A4 + A5 > 
(|23|) . Since unitarity and perturbativity constraints |A»| < 87r must hold (see for 
details [HI [TBI EH) verv large mass can arise only from very large negative 
m% 2 (EH)- The position on the (/ii , /U 2 ) plane of the actual state of the Universe, 
for the anticipated in the SM value of the Higgs mass Mh s < 200 GeV and 
Md± ^ 100 TeV, corresponds to /ii > 0,/i2 < 0, with large ratio of their 
absolute values > 10 5 - 10 6 , if Ai/A 2 ~ 1. 

In accordance with results of refs. [H QJjl [TJ, if A4 + A 5 > then the mixed 
phase M cannot exist (see e.g. ([2"510 while the charge breaking phase can. The 
charge breaking vacuum can be realized if in addition |i?3| < 1, where R3 — 
A3/1/A1A2. Simple analysis shows that the phase diagrams for this case are 
similar to those in Figs. I6.21 16.31 with the replacement R — > R3. Rays similar 
to the rays 41, 51, 52 and 54 give nothing new in comparison with the cases 
discussed in sect. [31 

The really new opportunity could appear for the ray similar to ray 53 in 
Fig. 16.31 (with rays going through the charge breaking vacuum). However, this 
opportunity is ruled out. Indeed, to realize it one needs C2 < and |c2 |/ci > 
l TO 22l/ TO ii ^ (10 5 H-10 6 ). The latter inequality contradicts the c\ > —c-i relation 
(|10p based on the positivity condition. It means that in our simple model 
Universe evolution to the current inert phase cannot pass through the charge 
breaking phase. 



8 Results and discussion 

Main results. The most important observation we made in this paper is as 
follows: if current state of the Universe is described by IDM, then during the 
thermal evolution the Universe can pass through various intermediate phases, 
different from the inert one. These possible intermediate phases contain no dark 
matter, which appears only at the relatively late stage of cooling down of the 
Universe. 

A complete set of possible ways of evolution of the Universe, including both 
EW symmetric (EWs) and EW non-symmetric (EW violating, EWv) initial 
states, can be summarized as follow^: 



'h 



EWs 
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rays 11, 21, 31, 41, 51 
A M A h ray 32 



A h 



ray 12 




(44) 



^Symbol I or II over arrows corresponds to the type of phase transition 
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We see that both a simple EWSB with a direct transition to the inert phase, 
as well as sequences of two or three transitions from EWs to the inert phase 
are possible. We found also that the current inert state of Universe can be 
obtained from both initial high temperature state with EW symmetry and from 
the initial state without this symmetry. 

As far the charged DM is concerned, still in principle allowed by the data, 
if heavy enough, it seems to be excluded in IDM. 

To find what scenario of evolution is realized in nature, one should measure 
all parameters of potential. The program how to measure these parameters at 
LHC and ILC is under preparation. 

Outlook. In contrast to the standard picture, these scenarios allow for the 
phase transition to the current inert phase at relatively low temperature, giving 
new starting point for calculation of a today's abundance of the neutral DM 
components of the Universe. 

In this paper we calculated thermal evolution of the Universe in the very 
high temperature approximation, i.e. for T 2 ^> |m|j|. The most interesting 
effects are expected at lower temperatures, where more precise calculations are 
necessary. The simplest expected modifications of the presented description are: 

1. Appearance of cubic terms like 3 T |17) . These terms are important near 
phase transition point, as they can transform some second-order phase 
transitions into the first-order transitions. 

2. The parameters become depend on temperature in more complicated way 
than that given by ©. Therefore, the rays, depicted thermal evolutions 
in Figs. 16.11 16.21 and 16.31 can become non-straight. The bending of these 
rays can be different in different points of our plots and at different . It 
can give possible spectrum of phase evolutions even reacher that discussed 
above. 

However, we expect that the general picture will not change too much. 

Possible extensions of model. The model we focused on in this paper con- 
tains two doublets. One can consider similar model with three doublets (partic- 
ular case of 3HDM) , as it is discussed in ref. [H] . It contains two standard Higgs 
doublets of 2HDM 4>s\ and 4>s2, coupled to fcrmions, and one Higgs doublet 
cf>D (having no coupling to fcrmions). The potential is invariant under S, D 
transformations: 

s s s 

{4>si, 4>S2} — > -{<f>si, 4>S2\, 4>d — > 4>d, SM — > SM; 

{4>si, 0S2> A {<f>si, S2 }, <\>d A -cf> D , SM A SM. (45) 

This model can incorporate all phenomena which appear in the standard 
2HDM, at the same time it contains DM particles as in the IDM. Thermal 
evolution of its parameters gives very diverse phase diagram which contain in 
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addition to the phases discussed above other phases, discussed in refs. [HO?] 
and their mixtures. However, even for this model our main conclusion about 
possible transformation of the Universe through the phase without DM holds. 

Reacher variants of both S- and D-sectors can be considered similarly. For 
example, one widely discussed model of this type (see e.g. [3D]) contains the 
same S'-sector as in our paper, but £>-sector consists of one doublet and one 
singlet scalars, non-interacting with fermions. This model has phenomenology 
similar to that discussed above but one hope to derive more or less natural 
values of couplings starting from SO(10) universality at the GUT scale. The 
biography of Universe in this model can be studied as above - obviously it will 
give more diverse phase story. 
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